Variations in the phenological patterns of a caddisfly inhabiting the same mountain massifs: Life‐history differences in different altitudinal zones

Abstract Organisms inhabiting mountainous regions can experience large vertical environmental changes, and show different ecological characteristics between altitudes, thus facilitating allopatric fragmentation even in geographically close populations. This study compared the life‐history patterns of a species of limnephilid caddisfly, Asynarchus sachalinensis, in several genetically differentiated populations between alpine and sub‐alpine zones in a temperate mountainous region. We showed that in the sub‐alpine populations, larval development started earlier with increasing water temperature in spring, and adult emergence was also earlier. The occurrence of adults was extremely low in mid‐summer, probably due to summer diapause, followed by a larger number of ovary‐developed females in autumn. On the other hand, in the alpine zone, increasing water temperature was delayed compared to the sub‐alpine zone, and larval development occurred from early to mid‐summer. Adult emergence and ovary‐developed individuals were concentrated in mid‐summer. Hence, summer diapause was not observed. These results indicated life‐history differences between genetically differentiated populations at different altitudes. As the timing of adult occurrence and ovarian developmental patterns differ between populations at different altitudes, it is possible that reproductive isolation is facilitated or maintained between populations.

. In addition, daylength, nutritional conditions and water temperature have been shown to be particularly important factors regulating life-history stages (postembryonic development) of aquatic insects that spend their larval stages in water (Bonacina et al., 2023;Füreder et al., 2005;Hayashi, 1990;Sweeney, 1984;Vannote & Sweeney, 1980;Ward & Stanford, 1982).Especially, in the case of lentic waters at high altitudes, the source of water in waterbodies is often dependent on snow cover, and the pattern of water temperature fluctuations varies considerably from place to place depending on the amount of snow cover, the location of the snow cover, and the timing of snow melt (Nebeker, 1971).Therefore, habitat differentiation is high, and life-history differences are correspondingly highly likely to occur between habitats (Butterfield, 1996;Goffová et al., 2015;McCulloch & Waters, 2018;Ward & Stanford, 1982).In addition, if there are differences in the reproductive timing of each population, reproductive isolation between populations drives speciation and/or intra-specific differentiation (Rundle & Nosil, 2005;Schluter, 2001).
For several cold-adapted species of limnephilid caddisflies, which inhabit lentic waters such as ponds and springs, it is known that they emerge with undeveloped ovaries and undergo diapause (i.e., "imaginal diapause") in caves and tree canopies until the temperature drops and/or water levels recover (Crichton, 1971;Denis, 1978;Novak & Sehnal, 1963;Salavert et al., 2008;Svensson, 1972).However, the occurrence and period of such imaginal diapause vary between species and populations (Crichton, 1971;Denis, 1978;Gíslason, 1978;Svensson, 1972).Factors contributing to differences in the occurrence and periods of imaginal diapause are explained by, for example, temperature and/or photoperiod (Alekseev & Starobogatov, 1996;Denis, 1978).One of the previous studies have shown experimentally that Limnephilus rhombicus, a species of limnephilid caddisflies, is prompted to induce and terminate "imaginal diapause" based on the daylength experienced during the 5th instar larval period (Denis, 1978).Differences in developmental timing caused by differences in water temperature may also occur, and along with differences in the daylength experienced by 5th instar larvae, these factors contribute to differentiation between populations and the occurrence of imaginal diapause (Bouchard & Ferrington, 2009;Finn et al., 2022;Hogue & Hawkins, 1991).If there is differentiation in the timing of the adult (flying stage), it results in differentiation in the timing of the reproductive stage.Therefore, gene flow between populations is restricted, and such gene flow and genetic differentiation are facilitated.In low-latitude temperate mountainous regions, it may result in differences in the life histories of the species inhabiting such areas because of large environmental differences between altitudinal zones (Laiolo & Obeso, 2017;Shah et al., 2017).However, there have been few studies comparing the life histories of the same species between altitudinal zones in such mountainous regions.
The target species of our study, Asynarchus sachalinensis Martinov 1914 (Figure 1), is a species of Limnephilidae with distributions in Sakhalin, Hokkaido and central to northern regions of Honshu.
Asynarchus sachalinensis inhabits cool temperature water bodies (e.g., ponds and springs) (Ito, 2008;Nozaki, 2002), and also high-altitude wetlands and caldera lakes of alpine zones (Nozaki, 2002;Suzuki et al., 2024).A previous life-history study of populations distributed in Hokkaido showed that larvae grow rapidly and synchronously in spring (April to May), emerge in early summer (in June), enter imaginal diapause in mid-summer (July to August), and lay eggs in autumn (September to October) (Ito, 2008).We showed in a previous study that this caddisfly has a univoltine lifecycle (Ito, 2008).In addition, our own previous study showed that A. sachalinensis collected from the entire distributional area including this study area (i.e., Mt.Norikura-dake, alpine and sub-alpine zones) consists of two genetically differentiated intra-specific clades based on mitochondrial DNA sequences (Suzuki et al., 2024).For these intra-specific clades, a relationship between the types of habitats and genetic differentiation was indicated (Suzuki et al., 2024).One of the genetic clades inhabits low lands while the other, inhabits high-altitudes area.Habitats at high altitudes are often isolated and discontinuous with nearby lowlands (e.g., basins).With such a pattern, the loss of population connectivity restricts chances for gene flow between populations (Knowles, 2001;Mikami et al., 2023;Pauls et al., 2006;Theissinger et al., 2013;Uscanga et al., 2021).In addition, even if populations are geographically close to each other, gene flow between populations is highly restricted, such that populations often show distinct genetic structures from each other (Mikami et al., 2023;Theissinger et al., 2013;Uscanga et al., 2021).
These considerations led us to study the relationship between habitat and life history in six populations (three populations inhabiting alpine zones and three populations inhabiting sub-alpine zones) of a cold-adapted species on Mt.Norikura-dake, one of the southernmost mountain ranges of the species' distribution.

| Study sites and a targeted species
Mt. 137.554 E,3026 m asl.), where the six water bodies studied are located, is a cluster of active volcanic calderas in the Chubu Mountainous Region of Honshu in the Japanese Archipelago (Shimizu, 1990) (Figure 2).The results of our previous study established that the water bodies of Mt.Norikura-dake are inhabited by large numbers of the target species Asynarchus sachalinensis at two different altitude zones (Suzuki et al., 2024).Previous studies indicated the forest limit to be at approximately 2500 m alt. on the east side of Mt.Norikura-dake (Miyajima & Takahashi, 2007;Takahashi et al., 2003).Of the two groups of study sites, one is located between approximately 1200 m and 1500 m asl.(SA1: Lake Ohmagari-ike, 1424 m asl.; SA2: Lake Chidori-ike, 1252 m asl.; SA3: Miyanohara wetland, 1227 m asl.), within the sub-alpine zone where deciduous broadleaved trees dominate the vegetation (Miyajima & Takahashi, 2007;Miyajima et al., 2007).The other study sites are at approximately 2700 m asl.(AL1: Lake Kiezuga-ike, 2737 m asl.; AL2: Lake Gono-ike, 2725 m asl.; AL3: Lake Tsuruga-ike, 2703 m asl.), located within the alpine zone where the alpine dwarf pine tree dominates (Miyajima & Takahashi, 2007;Miyajima et al., 2007).Dead broad-leaves were found on the bottom of water bodies in the sub-alpine zone and dead shrub cores (i.e., of alpine dwarf pine trees) were found on the bottoms of water bodies in the alpine zone.The water bodies of all study sites were lentic.The area (m 2 ) and circumference (m) of the water bodies were calculated based on 5 m mesh DEM maps from GSI (Geospatial Information Authority of Japan: https:// www.gsi.go.jp/ ).The information on the study sites is shown in Table 1, Figures 2c,d and S1.
Phylogenetic analyses conducted in our previous study using the mitochondrial DNA (mtDNA) combined sequences of the COI (648 bp), and the 16S rRNA (417 bp) regions showed that the respective populations of A. sachalinensis within these two altitudinal zones belong to distinct genetic clades, and that populations within the same altitudinal zone are genetically close (Figure S2: Suzuki et al., 2024).

| Sampling design for assessment of life-history differentiation between the alpine and sub-alpine zone populations
Weekly larval developmental stage surveys were conducted to compare the life histories of A. sachalinensis in those two different altitude zones of Mt.Norikura-dake (conducted surveys every week 32 times, from April 10 to November 17, 2022).In order to compare the larval period of A. sachalinensis inhabiting alpine and sub-alpine zones, larvae and pupae were collected quantitatively at each study site and on each sampling date.For collection of larvae, the same researcher selected five sampling points at each study site randomly and on each sampling date and collected by hand all larvae in quadrats (0.3 × 0.3 m).To identify larvae instar, image analysis software Image J 1.53 k (Schneider et al., 2012) was used to measure head width.Head width measurements and histogram construction using head width data were based on methodology used in a previous study of the life cycle of A. sachalinensis (Ito, 2008).Multiple histogram peaks were detected and the number of peaks was used to identify the instar stage number.
Adults were collected using light traps with 3 W ultraviolet LEDs (Figure S3) from sunset to sunrise at the same locations (two sampling points) at each study site.All adult specimens collected that had fallen onto the surface of the water in the light trap trays were counted while distinguishing between males and females.Adult emergence surveys were continued on two consecutive days at weekly intervals in each altitudinal zone until no adults were observed at all (the start observation date for the alpine zone was July 8, last observation was October 18; the start observation date for the sub-alpine zone was June 14, last observation date was November 17).However, one of the scheduled surveys in the alpine zone was not conducted due to bad weather conditions (October 6).To measure and compare environmental factors at each study site, EC (electrical conductivity), DO (dissolved oxygen), and pH were measured using an EC meter (Cyberscan CON 400, Nijkerk, Netherlands), a DO meter (MM-41DP, TOA DKK, Tokyo, Japan) and a pH meter (392R, AS ONE, Osaka, Japan).Surface water was collected at five points (from which larvae were collected) at each study site and on each TA B L E 1 Study site information list.sampling date, and each sample was measured for EC, DO, and pH on site.The averages of the measured data were used as the representative data recordings for each study site and sampling date.To assess the canopy openness rate, digital photographs were taken with a camera equipped with a fisheye lens (EX-FR200, CASIO, Tokyo, Japan) at six points around the ponds at each study site every month.Each photo was image-analyzed using CanopOn2 software (http:// taken aka-akio. org/ etc/ canop on2/ ), and then averaged.In the image analyses, 10% of the image margin was removed in order to eliminate the influence of artificial objects and/or background mountains.In adult emergence surveys, the air temperature was recorded hourly during the period from sunset on each survey day to sunrise on the following day (i.e., night-time air temperature) using data loggers (TRW-200WR, TOHO, Kanagawa, Japan).The water temperature at each study site was also recorded hourly using data loggers (MX2201, HOBO, Massachusetts, USA) during the same study period.A data logger measuring water temperature was settled at the bottom of each water body and hourly water temperature recordings were made throughout the study period, including the period of adult emergence.In a previous study, a relationship was indicated between daylength exposure during the 5th instar larval stage and the presence or absence of imaginal diapause.
When there were differences in larval development patterns of A.
sachalinensis between alpine and sub-alpine zones, there were differences in daylength exposure during the 5th instar larval stage.For this reason, daylength data was obtained from the National Astronomical Observatory of Japan (https:// eco.mtk.nao.ac.jp/ ) and correlated between the daylength and period of the 5th instar larval stage.

| Data analyses
Head width was measured on larvae collected on each sampling date and study site, and the dataset of all individuals (alpine zone: n = 324, sub-alpine zone: n = 911) was used to statistically test for differences in head width of larvae at each stage between alpine and sub-alpine zones.Data normality was tested using the Shapiro-Wilk normality test.The difference between larval head width of alpine and subalpine zones was analyzed by Mann-Whitney U-test.
Based on the number of days elapsed and the frequency of adult emergence observed using light traps with the starting date as "day 1," differences in the timing of adult emergence and ovarian development between the alpine and sub-alpine zones were statistically tested.
Data normality was tested using the Shapiro-Wilk normality test.The Based on the monthly average of each of the environmental factors (EC, DO, pH, water temperature, air temperature, and canopy openness) measured at each study site, principal component analysis (PCA) was conducted to investigate environmental differences between the alpine and sub-alpine zones (Table S1).
The Shapiro-Wilk normality test was performed using the function shapiro.test(), the non-parametric Mann-Whitney U test was performed using the function wilcox.test(),and PCA was performed using the function prcomp().These functions were implemented in statistical analysis software R4.2.2 (R Core Team, 2022).PCA plots were plotted using R package ggbiplot (https:// github.com/ vqv/ ggbiplot).

| Comparison of the environmental factors between the alpine and sub-alpine habitats of Asynarchus sachalinensis
The results of the PCA of the environmental factors showed that the first accounted for 64.0% of the variance and the second axis 15.5% (Table S2; Figure 3).The ellipsoid 95% probability confidence interval showed that habitats in the alpine and sub-alpine zones were largely differentiated environmentally (Figure 3).The factor loadings are shown in Table S2.The results showed that the environmental factors most highly contributing to variance of the first axis were canopy openness (−0.500),Water temperature (0.489), EC (0.488), and DO (0.450).The primary factor contributing to in second axis was pH (−0.813) (Table S2).During the adult emergence surveys, night air temperatures tended to be lower in the alpine zone than in the sub-alpine zone (Figure 4c).

| Larval development and the habitat environments of each altitudinal zone
No significant differences were indicated between the head widths of the five instars inhabiting the alpine and sub-alpine zones (Table 2, Figure S4).The patterns of developmental stage and adult emergence of A. sachalinensis did not indicate that different cohorts were present.Thus, A. sachalinensis was found to be univoltine in the populations of both altitudinal zones.Additionally, egg masses were observed under stones on land, indicating that A. sachalinensis overwinter as eggs or larvae in the egg masses during the winter.
Figure 5 shows the differences in the timing of the larval developmental stage and the pattern of water temperature fluctuations in the water bodies at each study site.Regarding the start of larval development, the populations inhabiting the alpine zone started later than the populations inhabiting the sub-alpine zone, and pupae were observed in the alpine zone more than 1 month later than those in the sub-alpine zone (Figure 5a,c).In particular, in the water bodies AL1 and AL2 in the alpine zone, there was a tendency for delayed larval development compared to AL3, where a significant rise in water temperature begins earlier in the season.The larvae in AL1 and AL2 had a shorter period in the fifth instar stage compared to AL3 (Figure 5a).
The date and daylength of first observation of fifth instar larvae at study sites was different in each altitudinal zone; SA1 and SA3 were April 24 and daylength was 13.46 h, while SA2 was May 1 and daylength was 13.62 h in the sub-alpine zone.In the alpine zone, AL1 was July 28 and daylength was 14.05 h, AL2 was July 13 and daylength was 14.38 h, and AL3 was June 21 and daylength was 14.55 h (Table S3, Figure 5).

| Adult emergence and ovarian development
Comparison of adult emergence patterns between the alpine and sub-alpine zones showed significant differences in the emergence patterns for all categories of all individuals, males, females, and females undeveloped reproductive organs (Table 3, Figure 4c).In the alpine zone, adults were observed intensively from the end of July to the middle of September.At each study site in the alpine zone, the dates when the largest number of adults were observed were August 24 for AL1 and AL2 and July 28 for AL3.In the sub-alpine zone, adults were observed from mid to late June, but almost no adults were observed in July and August, when air temperatures were higher.Adults were observed again from the middle of September to the middle of October, when air temperatures began to fall.The date when the largest number of adults was observed was September 27 at all study sites in the sub-alpine zone (Table S3, Figure 4).
The pattern of Stage B female emergence was compared between the alpine and sub-alpine zones.In the alpine zone, the emergence of Stage B females was observed approximately 20 days Females with spermatophores were observed continuously from mid-summer to early autumn in the alpine zone.In the sub-alpine zone, females with spermatophores were observed intermittently in early summer, mid-summer, and autumn (Figure 6b,d).

| The relationships between larval development and environment
In aquatic insects, it is well known that differences in larval development (e.g., ecological differences/life-history differences) occur due to differences in the altitude of habitats.In particular, habitat temperature is one of the most important environmental factors that alters life-history patterns (Bonacina et al., 2023;Füreder et al., 2005;Hayashi, 1990;Sweeney, 1984;Vannote & Sweeney, 1980;Ward & Stanford, 1982).
Several studies of limnephilid caddisflies have shown that the larval developmental patterns are determined by water temperature (Gíslason, 1978;Gotceitas & Clifford, 1983;Wissinger et al., 2003).The present survey showed that the season of appearance of early instar larvae of Asynarchus sachalinensis differed between the alpine and subalpine zones (alpine zone: early summer season, sub-alpine zone: spring season).The seasonal timing of larval development and the transition patterns also differed between the alpine and sub-alpine zones.
Populations of A. sachalinensis inhabit lentic waters in alpine zones; A. sachalinensis is generally known to inhabit lentic waters and springs (Ito, 2008).Data on water temperature fluctuations indicated that water temperatures increased more slowly in the alpine zone than in the sub-alpine zone.It was indicated that A. sachalinensis exhibited delayed development in the alpine zone waters due to the delayed increase in water temperature.
Additionally, the amount of snow cover varied between the study sites in the alpine zone, indicating that in AL1, where the snow cover was heaviest, the rise in water temperature was delayed and growth was delayed accordingly.However, fifth instar larvae were observed over relatively long periods in sub-alpine water bodies (SA1-3).In such water bodies, snow cover is relatively lower and snowmelt begins earlier in the season (Figure 5).As such, it was considered probable that factors other than water temperature may be involved in pupation initiation.
In particular, in the studies of the life-history of caddisflies, differences in daylength are known to be a factor in the induction of pupation and the initiation of larval diapause (Dixon & Wrona, 1992;Elliott, 1968;Malicky, 1981;Tsuruishi, 2003).Life-history adjustment is based on changes in daylength such that differences in developmental stages between individuals within a population are homogenized, and as a result, reproductive seasons are effectively synchronized (Wiggins, 1996).It is assumed that the daylength to which larvae are exposed affects the timing of emergence and the presence or absence of imaginal diapause, resulting in the synchronization of breeding timing within the population.In fact, our studies indicated that within each altitudinal zone, the emergence of adults was largely synchronous (Figure 4a,b).Therefore, it is indicated that the daylength experienced as larvae was related to their ecological characteristics (i.e., imaginal diapause) as adults.
In the sub-alpine zone, the first observation of the adults was in early summer (in June).However, few adults were observed during mid-summer (in July to August).Of these, a few adults observed shortly after emergence in early summer had not developed ovaries.Most of the individuals observed during autumn (in September to October) had fully developed ovaries (Figure 6c).These results Note: Females (Undeveloped reproductive organs): Developmental stage of ovaries A (Gower, 1963).
Females: Developmental stages of ovaries B to D (Gower, 1963).
Length of flight period (first date-last date) is recorded as the number of days specimens were collected in the light trap, beginning from the first day which is treated as "day 1." indicated the occurrence of imaginal diapause within sub-alpine populations of A. sachalinensis.To explain the decline in the number of adults in the sub-alpine zone in mid-summer, it was assumed that adults were not attracted to the lights due to imaginal diapause.
The presence of imaginal diapause in A. sachalinensis was already indicated by a previous study of populations in Hokkaido (Ito, 2008).
Our study showed similar ecological characteristics of populations inhabiting sub-alpine zones in Honshu.
On the other hand, in the alpine zone, the emergence of adults was concentrated in mid-summer (in July to August).Ovarian development of females was observed to be complete earlier than in the sub-alpine zone, about 30 days after the first adult emergence in the alpine zone.Therefore, imaginal diapause was assumed to be absent or shortened within alpine populations, in contrast to the sub-alpine populations.The lower temperatures in the alpine zone compared to the sub-alpine zone may have allowed A. sachalinensis, which is ancestrally adapted to cooler environments, to remain active during the summer season.A previous study comparing the life histories of Limnephilus affinis, which is a closely related species to A. sachalinensis, between populations at different latitudes and altitudes observed that at high latitudes and high altitudes in cooler environments, the adult flight period was concentrated within the summer season (Gíslason, 1978).
Previous studies on limnephilid caddisflies have shown a relationship between imaginal diapause and the daylength during the 5th instar larvae and/or adult stages (Denis, 1978;Novak & Sehnal, 1963).Asynarchus contumax, which is a closely related species to A. sachalinensis, is distributed at high latitudes of Palaearctic, matures during long-day conditions, and undergoes no imaginal diapause (Solem, 1983).For A. sachalinensis, whose life-history pattern differs between altitudinal zones, the larval development period likewise differs between the alpine and sub-alpine zones.Because larval development was at different periods between the alpine and sub-alpine zones, the daylength of exposure of 5th instar larvae differs even within the same mountain massif.Therefore, differences in the daylength experienced by 5th instar larvae are assumed to be one of the factors influencing the presence or absence of imaginal diapause.This would result in gene flow being highly restricted between these altitudinal zones, strongly supporting the maintenance of genetic differentiation between these zones.

| Genetic and ecological differentiation due to altitude
Habitats at high altitudes are often isolated and discontinuous with nearby lowlands (e.g., basins).With such a pattern, the loss of population connectivity restricts chances for gene flow between populations (Knowles, 2001;Mikami et al., 2023;Pauls et al., 2006;Theissinger et al., 2013;Uscanga et al., 2021).
Phylogenetic analyses based on sequencing of the mtDNA COI region conducted in previous studies showed differing genetic structures between alpine and subalpine Asynarchus sachalinensis populations.In particular, a single unique mountain massif-specific haplotype was detected in the alpine populations, which was not found among other populations (Suzuki et al., 2024).Therefore, such alpine populations are isolated and genetic differentiation has occurred (Suzuki et al., 2024).Not only geographical isolation, but also differences in life history between alpine and sub-alpine zones and differences in breeding season would lead to reproductive isolation.
In our study, the pattern of occurrence of females with spermatophores in their bursa is different between the alpine and sub-alpine zones.The females with spermatophores indicated that they had mated immediately after emergence within each altitudinal zone.
Thus, showing that the mating seasons correspondingly differ between the altitudinal zones.
No previous studies have provided data on the flight ability of A. sachalinensis.However, it has been observed that Chaetopterygopsis maclachlani (Limnephilidae), which inhabits the European Alps, has very low flight ability (Lehrian et al., 2010;Wagner, 1992).Phylogenetic analyses conducted based on sensitive molecular markers (e.g., genome-wide SNPs and/or microsatellite markers) will show that there is limited dispersal and gene flow between altitudinal zones.Therefore, the knowledge obtained from our other study will provide an understanding of the relationship between genetic and ecological differentiation between altitudinal zones.

F
I G U R E 1 (a) Asynarchus sachalinensis 5th instar larva crawling bottom of the water body (AL2).(b) A female adult (AL1).
E 2 Schematic maps of the study sites.(a, b) Mt.Norikura-dake (36.106N, 137.553E, 3026 m asl.)where the study sites were located is in the Chubu Mountainous Region in central Honshu.(c, d) A plain view map and a bird's-eye view map of Mt.Norikura-dake.Study sites AL1-3 are located in the alpine zone, and sites SA1-3 sites are in the sub-alpine zone of Mt.Norikura-dake.
The degree of ovarian development was classified into four categories-Stages A-D, based on Stage A: a stage with undeveloped reproductive organs (non-egg-maturation stage); Stage B: a stage in which yolk accumulation is insufficient (egg-maturation stage); Stage C: a stage in which yolk accumulation is sufficient, with fully mature eggs (egg-maturation stage); Stage D: a stage in which the ovary regresses after oviposition (post-oviposition stage).Furthermore, females with spermatophores in the bursa copulatrix were counted with reference to Ito (2020).
differences between adult emergence and ovarian development of the alpine and sub-alpine zones were analyzed using the Mann-Whitney U-test.For all collected adults at each study site, on each sampling date (alpine zone: n = 273, sub-alpine zone: n = 596), differences in adult emergence patterns between the alpine and sub-alpine zones were assessed for the categories of all individuals, males, females, and females undeveloped reproductive organs.Females with Stage A F I G U R E 3 Principal Component Analysis (PCA) of the study sites of Asynarchus sachalinensis based on the five environmental factors described in TablesS1 and S2.The points marked with red circles correspond to study sites located within the alpine zone.The points marked with blue squares correspond to study sites located within the sub-alpine zone.Ellipses indicate the 95% confidence intervals.ovarian development were treated as immature females and females with Stage B to D ovarian development were treated as mature females.However, the data about adult emergence collected on June 14 were not used in the analyses as it was only a pilot survey.Therefore, June 21 was set as the start date of the collection of light trap data used for the analyses.FI G U R E 4 (a, b) The number of adults of Asynarchus sachalinensis [Males, Females and Females (Undeveloped reproductive organs)] for each study site in the alpine and sub-alpine zones.(c) The trends of the average night air temperature at each study site and their respective daylength at Mt. Norikura-dake.

6
Number of females of Asynarchus sachalinensis observed in the alpine and sub-alpine zones.(a, c): The number of females observed in each altitude zone, classified according to their stage of ovarian development.The colors correspond to the color scheme of their ovarian development classification stages in (e).(b, d): The number of females observed with spermatophores within each altitude zone.(e): Typical stage ontogenetic patterns when ovarian development classified into stages A-D based on Gower (1967).